arXiv:astro-ph/0512651v3 21 Aug 2008 


Inhomogeneity-Induced Cosmic Acceleration in a Dust Universe 


rQ Je-An Gujil 


Chia-Hsun Chuangp Je-An GujU and W-Y. P. Hwan^ 

Department of Physics, National Taiwan University, Taipei 10617, Taiwan, R.O.C. 

Abstract 
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I. INTRODUCTION 


In 1998, via the distance measurements of type la supernovae, it was discovered that the 


expansion of the universe is accelerating 


of 


The accelerating e^^ansion of the present 


universe was reinforced recently by the updated supernova data (3|, Ij, y, la, l7| and the WMAP 
measurement Q| of cosmic microwave background (CMB). Regarding the cosmic evolution, it 
is the common consensus that normal matter (such as protons, neutrons, electrons, etc.) can 
only provide attractive gravity and therefore should always slow down the cosmic expansion, 
i.e.. 

Normal Matter ^ Attractive Gravity ^ Deceleration. (1) 

Thus, to explain this surprising, mysterious phenomenon of the accelerating expansion, 
many people rely on exotic energy sources, as generally called “dark energy”, which provide 
signihcant negative pressure and accordingly anti-gravity (repulsive gravity). 

The above conclusion about the existence of the cosmic acceleration and the necessity 
of introducing dark energy is based on a simplihed cosmological model, the Friedmann- 
Lemaitre-Robertson-Walker (FLRW) model, and indeed could be model-dependent. In the 
FLRW model the universe is assumed to be homogeneous and isotropic (i.e. the cosmological 
principle) and accordingly the Robertson-Walker (RW) metric is invoked in the Einstein 
equations that describe the evolution of the universe. Meanwhile the energy-momentum 
tensor in the right-hand side of the Einstein equations is regarded to truly reflect the real 
energy distribution (averaged in space) of the universe. Nevertheless, so far there is no 
convincing proof validating this simplihcation. 

Apparently, our universe at present is not homogeneous at small scales. The cosmological 
principle is roughly realized only at very large scales. To take advantage of the cosmological 
principle and invoke the RW metric, it is necessary to perform spatial averaging over large 
scales, along with which the form of the Einstein equations in general should change because 
of the non-linearity of the Einstein equations 9(]. That is, (i) when invoking the RW metric 
and the real energy sources of our universe in the Einstein equations for describing the cosmic 
evolution, the left-hand side (the geometry part) of the Einstein equations should be modi- 
hed, or, from another point of view, (ii) if one insists to use the Einstein tensor corresponding 
to the RW metric in the left-hand side of the Einstein equations (or, equivalently, moving the 
above mentioned modihcation in the geometry part to the right-hand matter-energy part). 
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there should appear new, effective energy sources coming from geometry (which is certainly 
“dark”) and consequently the energy-momentum tensor in the right-hand side does not truly 
correspond to the real energy distribution of our universe. Thus, generally speaking, it is 
doubtful to employ the Einstein equations to describe long-time, large-scale phenomena, 
such as the evolution of our universe, while the spatially averaged energy-momentum tensor 
and the spatially averaged metric tensor are used therein. (For more discussions about the 
validity and the problems of the FLRW cosmology, see [f^.) 

nstead of invo king exotic energy sources or unconfirmed physics, it is suggested Q, 

3, H, H, H, HQ, Q, Q, 


ll3|. 


22l | that the cosmic acceleration might originate from the 


violation of the cosmological principle, homogeneity and isotropy, i.e., the acceleration might 


be induced by the inhomogeneities of the universe. The possible change of the dece 
parameter in an inhomogeneous universe has been pointed out in Refs. l23| and 


eration 


24. In 


particular, it has been shown that the luminosity distance-redshift relation indicated by 


the supernova data can 


De reproduced in an inhomo gen eous cosmological model without 


introducing dark energy [25|, l2a, l27|. 


28|, 


291, lad, laii, i32|. 


aai. 


m- 


The current situation of our knowledge (what we know and what we do not know well) 
about observations, cosmic acceleration and inhomogeneities is as follows. 


— Known — 

• Based on the FLRW cosmology, the current observational results indicate the existence of 
the cosmic acceleration. (That is, by invoking the FLRW cosmological model to interpret 
the observational data, one would conclude that the expansion of the universe is accelerating 
in the recent epoch.) 

— Doubts — 

(1) Is the FLRW cosmology a good approximation? 

(2) Do the current observational results indicate the existence of cosmic acceleration for 
the real universe with complicated energy distribution? (Even the definition of accelerating 
expansion is an issue for our complicated universe. This issue will be discussed in Sec. ttH) 

(3) Can the inhomogeneities of our universe explain the observational results? 

(4) Can the inhomogeneities of our universe generate accelerating expansion? 

These four doubts are far from being fully answered. This is due to the difficulties from 
the complexity of the real energy distribution and the non-linearity of the Einstein equations. 
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Because of these difficulties, instead of dealing with the full Einstein equations describing 
the real complicated universe, usually people employ the following two approaches to sketch 
the possible answers to the above doubts. 

(I) Taking perturbative approach with a simple background (such as the homogeneous and 
isotropic RW background). 

(II) Invoking exact solutions of the Einstein eguations describing an inhomogeneous universe. 


In addition to these two approaches, there are also non-perturbative studies without invoking 


exact solutions, e.g., see Refs. l35l. 


36 


37 


Via Approach (I) with a homogenous and isotropic background, the positive answer to 


Doubt (1) is supportec 
supported in Refs. [3^, 


in Refs. 38 


39 


40,141. I 


and 


39|, and the negative answer to Doubts (3) and (4) 


Nevertheless, the reliability of the perturbative approach 


in (I) is doubtful for investigating the late-time cosmic evolution, for which no 


distribution, but also curvature, may not be described with perturbations 


more, the arguments of Refs. 


38 


39, 


40 


41 




have been countered in Refs. 


21 


only energy 


421. Further- 


313 , 


43. 


The drawback of Approach (II) is that the exact solution invoked may be very different 
from the real situation of our universe. From some angle, instead of Doubts (3) and (4), 
Approach (II) is to answer a more general question: 

• Can inhomogeneities explain observational results and generate accelerating expansion? 
Or, more precisely, 

(3') Does there exist a universe (maybe very different from ours) in which the inhomogeneities 
can explain the observational results? 

(4') Does there exist a universe in which the inhomogeneities can drive the expansion to 
accelerate? 


Via Ap proach (II), for Doubt (3') it has been shown in Refs. 


33 


34 


25 


M 


U 


28 


29 


, |3Q|, 131|, 132, 


that the sup ernova data can be explained by invoking the Lemaitre-Tolman-Bondi 


solution 


M, 


45l. l46l | that describes a spherically symmetric (but inhomogeneous) dust fluid. 


In the present work, we take Approach (II) for investigating the general possibility of 
generating accelerating expansion via inhomogeneities, i.e. Doubt (4'). Against the common 
intuition in Eq. ([T]) we hnd examples of accelerating expansion in the case of a spherically 
symmetric dust fluid described by the Lemaitre-Tolman-Bondi (LTB) solution, thereby giv¬ 
ing support to the positive answer to Doubt (4'). In addition to our examples, Kai et ah [T7| 
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also found acceleration examples based on the LTB solution, where there exists singularity 
around the center of the spherically symmetric system during the accelerating epoch. In 
contrast, in our examples the system is smooth everywhere and no singularity is involved. 


(For more comparison, see Sec. 0) The possibility of the acce’ 
model was also pointed out by Paranjape and Singh in Ref. 


erating expansion in the LTB 


48l. in which numerical mod¬ 


els exhibiting acceleration were constructed by an approximation where the contribution of 
3-curvature dominates over the matter density. 

There a 

H Q, 

smoothly and the effect of the junction between these regions are the essential issues yet 
to be seriously explored. We particularly note that it needs much caution to connect two 
(or more) regions. In many cases the effect from the junction between separate regions is 


so existed acceleration examples for a system consisting of two or more regions 
1^ . 1^ . 1^ . In these acceleration examples, connecting the separate regions 


signihcant anc 
Singh in Ref. 


lould not be ignored. A similar doubt was also raised by Paranjape and 


52l . As a demonstration of how things may go wrong when the connection 
or the junction is not appropriately taken care, in Appendix [O we investigate in detail the 


acceleration examples studied by Nambu and Tanimoto in Ref. l49|, and show that actually 


there should be no acceleration in those cases when we properly connect the separate regions 
and seriously take the effect of the junction into account. 

The acceleration examples we hnd could be far away from the real situation of our uni¬ 
verse. We are not proposing to employ these mathematical examples or the LTB solution 
to describe our universe. These acceleration examples, which provide positive answer to 
Doubt (4'), are to demonstrate how inhomogeneities can drive the expansion to accelerate, 
thereby showing how our intuition about the interplay of gravity and the cosmic evolution 
may go wrong [i.e. against the common intuition in Eq. ([1])]. Accordingly, the effects of 
inhomogeneities on the evolution of the universe should be carefully restudied. 

This paper is organized as follows. In Sec. m a tricky issue of the dehnition of acceler¬ 
ation is discussed and two dehnitions to be utilized for searching acceleration examples are 
introduced. In Sec. IIIH the LTB solution is described. In Sec. HV] and Sec. |Vl we present the 
examples of the accelerating expansion corresponding to these two dehnitions of acceleration, 
respectively. A summary and discussions follow in Sec. IVTl 

Throughout the present paper we will use the units where c = SttG = 1. We note that in 
this unit system there is still one unit unspecihed, and, as a result, the value of one of the 
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dimensionful quantities can be arbitrarily set. For example, in the acceleration examples we 
will present, the physical size of the spherical region under consideration can be any length 
long (such as 1 fermi, 1 cm, 1 pc, 1 Mpc, Hubble length Hq^, etc.). Once the value of one of 
the dimensionful quantities (which are independent of c and G) is settled, all the units for 
the dimensionful quantities studied in the present work are specihed. 


II. DEFINITIONS OF ACCELERATION 


In cosmology, expansion and acceleration of the universe are a subject associated with 
the evolution of the space (in size), and should have nothing to do with the particle motion 
relative to the space. How to dehne the speed and the acceleration purely corresponding to 
the evolution of the space, meanwhile avoiding confusion and interference from the particle 
motion relative to the space? This is a tricky issue for a real universe with complicated 
energy distribution [lOt]. 

To characterize the evolution status of the space, usually one needs to invoke two quan¬ 
tities: a length (distance) quantity L and a time quantity t, with which L and L (where the 
over-dot denotes the time derivative) present the speed and the acceleration of the expansion, 
respectively. The tricky issue mentioned above then corresponds to the problem of making 
the choice of these two quantities. With different choices one has different dehnitions. It 
will be puzzling if the description of the space evolution (expansion and acceleration) status 
is not universal but depends on the dehnition one invokes. 

If it is inevitable to make the choice which the description of the evolution status of 
the universe is based on and maybe sensitively depends on, what is the reasonable physical 
choice? The very guiding principle for making the choice is that we are quantifying simply 
the evolution of the space and we need to make sure that the dehnition we invoke does 
not involve particle motion relative to the space nor something fake stemming from an 
inappropriate frame choice. 

In this section we will introduce two dehnitions, as to be called “/me acceleration” and 
'^domain acceleration ”, respectively involving two length quantities: (i) the distance between 
two points in space and (ii) the size of a domain in space. For these two kinds of lengths 
the above tricky issue relates to the choice of two points for the line acceleration and the 
choice of the domain (i.e. its boundary) for the domain acceleration. With diherent choices 
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(corresponding to different definitions) one might obtain very different conclusions about 
the acceleration/deceleration status. 

In general, with a certain choice of frame (coordinate system), the two spatial points 
and the boundary mentioned above can be hxed in the spatial coordinate space. With an 
improper choice of these two points and the boundary, or, equivalently, with an improper 
choice of frame (in which the spatial coordinates of the two points and the boundary are 
hxed), the change of the length quantity with time may be incapable of representing the 
true evolution of the space, which could be mixed up with the evolution of the frame. The 
contribution from the frame evolution to the acceleration of the chosen length will be called 
“frame acceleration” in the present paper. 

To construct suitable dehnitions of acceleration for these two kinds of lengths, for simplic¬ 
ity we consider a universe consisting of freely-moving particles (i.e., moving along geodesics), 
among which there is no interaction other than gravity. In this case, for the line acceleration 
a simple reasonable length quantity to choose is the distance between two (freely-moving) 
particles. In contrast, an apparently improper dehnition which could lead to fake frame 
acceleration is to invoke the distance between two points which move “outward” relative the 
particles in between, i.e., accordingly, there are more and more particles between these two 
chosen spatial points. 

As to the domain acceleration, a reasonable choice is the size of a spatial domain in which 
the number of particles is constant in time. In contrast, an apparently improper dehnition 
which could lead to fake frame acceleration is to invoke a domain whose boundary moves 
outward relative the particles therein, i.e., accordingly, there are more and more particles 
within this domain. 

For implementing these two acceleration dehnitions, meanwhile satisfying the above men¬ 
tioned requirements for avoiding the confusion from particle motion and fake frame acceler¬ 
ation, it is particularly benehcial to use the synchronous gauge. Note that in a dust universe 
the synchronous gauge can be chosen if and only if the vorticity vanishes. In the synchronous 
gauge the line element is as follows. 

= —dt^ + hij{:K, t)dx^dx^ , (2) 

where t is the cosmic time. In this gauge, the cosmic time t is simple and universal to choose 
in dehning acceleration. Regarding the length quantity, in this gauge the above mentioned 
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requirements are easy to meet because a point fixed in the spatial coordinate space is a 
geodesic. In particular, for a dust fluid with the energy-momentum tensor 

= pu>^u '^, = (1, 0, 0, 0), (3) 

the distance between two fixed points and the size of a domain with its boundary fixed in this 
coordinate space are simple and direct choices satisfying the requirements. In the following 
we will focus on this simple case of a dust fluid described in the synchronous gauge, and 
introduce two definitions of acceleration involving these two length quantities respectively. 


A. Line Acceleration 


Regarding two points in space with the proper distance L{t) between them at time f, 
it is reasonable to use L{t) and L{t) to characterize the expansion/collapse status and 
the acceleration/deceleration status of the space in between. The expansion rate and the 
deceleration parameter for the proper distance L{t) are defined in the usual way, as follows: 


Hl = IjL, 

LIL 'll 

^ ~l!l ~ 


(4) 

(5) 


The condition {Hl > 0, < 0} corresponds to the accelerating expansion of the proper 

distance between these two points in space, which is dubbed “line acceleration” in the present 
paper. We have found examples of the line acceleration in a dust universe, as to be shown 
in Sec. HVl 


B. Domain Acceleration 


Another definition of acceleration, as dubbed “domain acceleration” in the present paper. 


11 


21 


35|,l3^,l37, 


has been widely used in the literature 
a spatial domain D with a finite volume 

Vn = [ Vhd^x, 


42 


49 


50 


51 


52 


53 


54|. It is for 


( 6 ) 


>D 


where h is the determinant of the spatial metric tensor hij. Invoking the length scale of the 
domain, 

Ld = W. (7) 
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one can define the expansion rate and the deceleration parameter of the domain in the usual 
way, as follows: 


Hd = Ld/Ld , 

_ Ld/Ld 
= - 


LnL 


D-^D 


f 2 


( 8 ) 

(9) 


The condition {Hd > 0, gD < 0} corresponds to the accelerating expansion of the domain 
(in size), i.e., domain acceleration. 


As shown in 


55 


56 


571 . l58| |. for an inhnitesimal domain in a dust universe without 


vorticity the deceleration parameter qo is always positive, i.e., corresponding to local de¬ 
celeration. Nevertheless, the non-local deceleration/acceleration status of a domain with a 
nonzero hnite volume (in particular, the observational universe of the Hubble size) may be 
very different So far there is no no-go theorem excluding the possibility of negative qd- 
On the contrary, we have found examples of the domain acceleration in a dust universe, as 
to be shown in Sec. El 

We note that, in the special case of a homogeneous and isotropic universe described by the 
RW metric with the scale factor a{t), the expansion rates and the deceleration parameters 
dehned above are the same as those in the standard cosmology: 


H = a/a, 
a/a 

IP 


q = 


aa 


( 10 ) 

( 11 ) 


In the next section we will introduce the LTB solution, based on which we hnd examples 
of acceleration. (Reminder: the units where c = 8nG = 1 will be employed.) 


III. LEMAITRE-TOLMAN-BONDI (LTB) SOLUTION 


The LTB solution 


44 


45l. l46[ | is an exact solution of the Einstein equations for a spheri¬ 


cally symmetric dust fluid. The metric is given by 


,2 _ 1.2 , (Rir) dr" D2jn2 


ds = —dt + 


+ R^dn^ , 


( 12 ) 


l + 2E(r) 

where R is a function of the time coordinate t and the radial coordinate r, E{r) is an 
arbitrary function of r, and R,r denotes the partial derivative of R with respect to r. With 
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this metric the Einstein equations can be reduced to two equations: 



2 


p(t,r) 


2E{r) ^ 2M(r) 
2M'(r) 


( 13 ) 


(14) 


where M(r) is an arbitrary function of r and the over-dot denotes the partial derivative 


with respect to t. The solution of Eq. flT^ can be written parametrically by using a variable 
7] = J dt/R, as follows. 



where tb{r) is an arbitrary function of r. Summarily, there are one dynamical field, R{t,r), 
and three arbitrary functions, E{r), M{r) and tf,{r). For a given set of the three functions 
{E'(r), M(r), 4(r)}, we have a solution R(t,r) specified by Eqs. ffTSj) and ffT6|) . 

Regarding the behavior of the dynamical held and the three functions introduced above, 
it is reasonable to consider the requirements that there is no hole and no singularity in space 
and the energy density is non-negative and hnite. 

(1) For no hole at the center (i.e. the area of the spherical surface at r = ro goes to zero 
when ro goes to zero), R(t, r = 0) = 0. 

(2) For no singularity, we consider R,r {t, r) ^ 0. 

(3) For a non-negative and hnite energy density, 0 < p(t,r) < oo. 

According to these requirements, the dynamical held and the functions in the LTB solution 
should satisfy the following restrictions: 


( 17 ) 

( 18 ) 


r = 0 : R, R, M, E = 0 , 
r > 0 : RR,r > 0 , MM' > 0 . 


Without losing generality, in our study and in the remaining of the present paper we choose 


R, R,r , M, M' > 0 . 


( 19 ) 


From the relation 



( 20 ) 
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as derived from the Einstein eqnation ffT3|) . the sign choice in Eq. implies 


R<0. 


( 21 ) 


For more details abont the restrictions on the LTB solntion and the corresponding featnres 
of the dynamical field and the fnnctions presented above, see Appendix Note that the 
necessary and snfficient conditions of no shell crossing in a period of time are presented by 
Hellaby and Lake in Ref. 1^. In the present paper we consider the absence of shell crossing 
at some time, and accordingly the above-mentioned conditions are necessary conditions. 

By introdncing the following variables 


a{t,r) = 


R{t,r) 


k{r) = 


2E{r) 


Po{r) = 


6M{r) 


( 22 ) 


the line element in Eq. fll2l) and the Einstein eqnations flT^ and flT^ can be rewritten in a 
form more similar to that of the RW metric: 

a,rr\2 


ds = —dt -I- a 


1 + 


1 — k{r)r‘^ 


+ r^dVti, 


(23) 



p(t,r) 


k{r) Po(^) 

a? 3a3 

3a‘^r^{ar)^r 


The solntion in Eqs. flT^ and flTBD becomes 


a(r),r) 




Pojr) 

6k{r) 


1 — cos 




+ h(r), 


(24) 

(26) 


(26) 

(27) 


where ^ r/r — / dt/a . 

Now the dynamical field describing the evolntion of the space is a{t, r) and the three 
fnnctions to be given for specifying a solntion a{t, r) are /c(r), pair) and hir) [corresponding 
to E(r), M{r) and R{r), respectively]. Apparently, when all the fnnctions a{t,r), k{r), 
Po{r) and h^r) have no dependence on the radial coordinate r, we come back to the RW 
metric from Eq. fl23l) . the Friedmann eqnation from Eq. fl2Tl) and the formnla realizing stress 
energy conservation from Eq. ([25]). From the comparison with the RW metric, we can get a 
rongh pictnre of the LTB metric: a{t, r) plays the role of a spatially varying, time-dependent 
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scale factor describing the evolution of the space, k{r) corresponds to the spatial curvature, 
pair) relates to the physical energy density p(t,r), and tfe(ri) can be regarded as the initial 
time of the big bang at r = ri, i.e., the space at r = ri starts to expand from singularity 
[a{tb,ri) = 0] at the time t = tb(ri). 

In search of examples of accelerating expansion in a dust universe described by the LTB 
solution, we tried a variety of LTB solutions corresponding to different choices of the func¬ 
tions {k{r), pQ{r),tb{r)}, and eventually found examples among these tedious trials, as to 
be studied in the following two sections. We note that there is redundancy in the choices 
of the functions {E{r), M{r),tb{r)} or {k{r), po{r),tb{r)}. For example, for a monotonically 
increasing function M(r) = po{r)r^, one can choose, without losing generality. 

Pair) = constant, (28) 

which is the choice invoked in our search for the acceleration examples. 


IV. LINE ACCELERATION INDUCED BY INHOMOGENEITIES 


Regarding the demonstration of how the acceleration can be induced by inhomogeneity, 
naively we have better chance with larger inhomogeneity which may induce more signihcant 
acceleration. In the LTB solution, due to the spherical symmetry, the inhomogeneity lies 
along the radial direction, while there is no inhomogeneity and also no acceleration [as 
implied by Eq. (El])] in the angular directions. Accordingly the possible line acceleration 
induced by inhomogeneity must involve the length component in the radial direction. For 
simplicity, we focus on the proper distance between the origin (r = 0) and the point at 
r = at time t: 

fTL 

(29) 
Jo 

where 


Qrr — 


Jo 

{a + a,rr) 


1 + 2 E{r) 1 — k{r)r‘^ 

The deceleration parameter corresponding to this radial proper distance is 


Qr = 


TJ—jf' 


(30) 


(31) 


The sign of the deceleration parameter is determined by the sign of or, more precisely, 
the integral of from the origin r = 0 to r = i.e., d^^/g^dr. 
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In the LTB solution, inhomogeneity can be introduced by choosing inhomogeneous func¬ 
tions for k{r), po{r), and tb{r). In this section, for simplicity we introduce inhomogeneity 

through only k{r) by employing the following function: 

_ {hu + l){r/rkT>^ ^ 

1 -k (rjrkYf^ 

(32) 

while choosing 


Poif) = constant. 

(33) 

h{r) = 0. 

(34) 


The behavior of the function k{r) in Eq. fl32D is illustrated in Fig. [H For a large power 
this k{r) function mimics a step function with violent change around r = Vk, accompanying 
which large inhomogeneity is introduced. We note that when k{r) = 0 = E{r) [with 
arbitrary poir) and tb{r)] there is no line acceleration (see Appendix [B]), and therefore the 
inhomogeneity in k{r) seems to play an essential role in generating accelerating expansion. 

rk 

1 

_ 0 

-hk 

rk 

r 



-hi 


FIG. 1: The plot of the function k{r) invoked in the search for acceleration examples. 

With the above choice of the three arbitrary functions in the LTB solution, we have six 
free parameters to tune: 

{t,rL,po,rk,nk,hk). (35) 

In search of examples of the line acceleration in the radial direction, we surveyed this six- 
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dimensional parameter space and did find examples eventually^ Table [I] presents one of the 
examples with significant acceleration, —— 0.834 ~ 0{1). We note that, as suggested by 
the observational data with the analyses based on the FLRW cosmology, the deceleration 
parameter of our present universe is of the order unity and is negative in sign. 


TABLE I: One example of the line acceleration in the radial direction. 


t 

TL 

ho 

Tk 

nk 

hk 

Lip 

Ijp 

Ijp 

Qr 

1 

1 

1 

0.7 

20 

1 

1.29 

1.16 

0.868 

-0.834 


The further details of the example in Tabled are illustrated in Fig. [2l which presents the 
energy density distribution and the acceleration/deceleration status of the radial and the 
angular line elements. As shown in this figure, the spherically symmetric dust fluid consists 
of three regions: two roughly homogeneous regions — the inner over-density region with 
positive k{r) and smaller a{t,r) and the outer under-density region with negative k{r) and 
larger a{t,r) — and one transition or junction region, where the inhomogeneity locates, of 
those two homogeneous regions.^ In the inhomogeneous region with significantly changing 
energy density we have acceleration in the radial line elements, i.e., > 0, while in 

the homogeneous regions with smoothly distributed energy density we have deceleration. 
This result strongly supports the suggestion that inhomogeneity can induce accelerating 
expansion. In addition, as expected there is always no acceleration in the angular directions 
along which everything is uniformly distributed.^ 

For a demonstration of the scales of the size and other quantities of this system, in Table 
II we fix the only one unspecified unit by using the length unit: O.lMpc, IMpc, lOMpc 
and lOOMpc, and present the values of several dimensionful quantities, respectively. Among 
these quantities, t corresponds to the time under consideration, Lr the size of the system, 
p{r = 0) roughly the energy density of the inner region and p{r = tl) roughly that of the 
outer region. As shown in this table, regarding the same example in Table [H when the size 
of the system increases by one order of magnitude, the time increases by one order and the 

^ For 1 — fc(r)r^ to be positive, we consider two different sufficient conditions, {nu > 0, /ife > 1, Vk < 1} and 
{hk > —l,r< tl < 1}, and restrict our search in the cases satisfying one of them. 

^ Note that in the inner region around the origin r = 0 the energy density distribution is fiat and therefore 
there is no singularity and, moreover, no cusp behavior (or weak singularity) in this acceleration example. 
^ For a proof, see Appendix El The result is in Eq. (lAlip . 
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energy density decreases by two orders of magnitude. For this example to be consistent 
with the situation of our present universe, the time t should be ~ years and the energy 
density of the outer region should be similar to the average energy density of the present 
universe, ~ 10“^®g/cm^. One can see that these two conditions cannot be simultaneously 
satished in this example, and therefore this example by itself alone cannot describe the 
present universe. 

TABLE II: Corresponding to different length units, the values of several dimensionful quantities 
for the acceleration example in Table |T] are presented. 


Length unit 

Lr(Mpc) 

t(year) 

p{r = 0) (g/cm^) 

p{r = ri) (g/cm^) 

0.1 Mpc 

0.129 

3.26E5 

1.19E-29 

6.80E-33 

1 Mpc 

1.29 

3.26E6 

1.19E-31 

6.80E-35 

10 Mpc 

12.9 

3.26E7 

1.19E-33 

6.80E-37 

100 Mpc 

129 

3.26E8 

1.19E-35 

6.80E-39 


To study the dependence of the deceleration parameter on the six parameters in Eq. 
(I35|) . we use the example in Table [I] as a reference and tune one of the six parameters at one 
time, while keeping the other hve unchanged (i.e., with the values in Table [11). The results are 
shown in Fig. [3l The plot of qr versus Uk shows that we have larger acceleration for larger Uk 
that corresponds to larger inhomogeneity. This result again supports the possibility of the 
inhomogeneity-induced acceleration, in which, naively, inducing larger acceleration requires 
larger inhomogeneity. In addition, we have acceleration even for a moderate power (e.g., 
rifc = 5), that is, for the purpose of generating the line acceleration we do not need very 
large inhomogeneity. 
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FIG. 2: The plots of the physical energy density p and the quantities, and which 

characterize the local acceleration/deceleration status in the radial and the angular direction, 
respectively, for the example in Table HI We note that when r = ri = 1, = —0.066 and 

d‘iV^e = -0.095. 
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FIG. 3: Illustration of the dependence of the deceleration parameter qr on the parameters 
PQ,rk,nk, hf^), using the example in Tabled] as a reference (denoted by the large dot). 

V. DOMAIN ACCELERATION INDUCED BY INHOMOGENEITIES 


In this section we investigate the domain acceleration for a spherical domain, 0 < r < r^, 
with the volume 


rro 


Vd = dvr 




(■ro 


: dr = dvr 


a‘^r‘^{a + a,r r) 


dr. 


(36) 


/o yi + 2E(r) to .yi - k{r)r^ 

The length invoked for the domain acceleration is that in Eq. ([7]) : Ld = via which the 
deceleration parameter qo of the domain is defined in Eq. dd]): Qd = —LdLd/■ We note 


that in t 
and Ref. 


le case where k{r) = 0 = E{r) there is no domain acceleration (see Appendix | 


m- 
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For the domain acceleration Kai et al. ^ found examples based on the LTB solution, 
where constant poi^), trivial tb{r) and the following k{r) function are invoked. 


k{r) = < 


ko 

kp 

2r2 

kp 

2r2 

kp 

2r2 

0 


rt 


r 2'|2 

^ + rf + rl 


1 ' 2 
(r? + ri) 

(r? + ri) 



- 1 


for 0 < r < ri, 
for ri < r < r 2 , 

for r 2 < r < ra, 

for ra < r < rb, 

for Tb < r, 


( 37 ) 


where 0 < ri < r 2 < ra < rb and kp is a constant. In these examples 47(], the acceleration 
involves the existence of a singularity around the origin. In contrast, we hnd the domain 
acceleration examples without singularity. The difference in the choice of the three functions 
which specify the LTB solution is that in our examples non-trivial tb{r) is invoked, as going 
to be presented in the following. Our examples indicate that the existence of a singularity 
is not necessary for generating the domain acceleration. 

In search of examples of the domain acceleration without singularity, we first follow 
the same procedures used in the previous section for the line acceleration and choose the 
functions, fc(r), Po{r), and 4(r), as those in Eqs. fl52]) - fl5T|) . In this case we surveyed the 
six-dimensional parameter space (t, r r, P n, rk, rik, hk) and found no domain acceleration. 

Contradicting our result, in Ref. it was claimed that the example of the domain 
acceleration was found with constant pp, trivial tb{r) and step-function-like k{r) [i.e., with 
inhnitely large Uk in Eq. fl32l) ]. There is a mistake in the calculations of the volume Vd in 
Ref. where the authors ignored the volume at the transition point r = that is actually 
nonzero and should not be ignored (even though r = corresponds to a “2D surface” in 
the coordinate space). After taking the volume at r = back into account, we found no 
domain acceleration. (For more details, see Appendix [Cl 1 

Since no domain acceleration (without singularity) was found in our search with trivial 
tb{r), we then consider a non-trivial function for tb{r) : 

Kb{r/rtY^ 


hir) = 


1 -F {r/rtY^ ’ 


(38) 


18 











while invoking the same functions for k{r) and poir) in Eqs. (I3^ and fl33|) : 


k{r) 

Po{r) 


{hk + l)(r/rfc)^^ 

1 + {rlrk)'^k 

constant. 


The behavior of the function tb{r) in Eq. (138|1 is illustrated in Fig. 01 


( 39 ) 

( 40 ) 



r 


FIG. 4: The plot of the function invoked in the search for the domain acceleration. 

With the above choice of the three functions involved in the LTB solution, we have nine 
free parameters to tune: 

{t, td, Po, Tk, Uk, hk, n, nt, htb ). (41) 

We surveyed this nine-dimensional parameter space and did hnd examples of the domain 
acceleration eventually. In Table lllll we present three examples with signihcantly different 
magnitude in acceleration (i.e., regarding the value of qr)). 

The further details of these three examples are respectively illustrated in Figs. [5l [6l and 
[71 which present the energy density distribution and the acceleration/deceleration status of 
the radial and the angular line elements. Similar to the example of the line acceleration 
in the previous section, in all these three examples the spherically symmetric dust fluid 
consists of two roughly homogeneous regions and one inhomogeneous transition/junction 
region.^ Regarding the line elements, we get acceleration in the radial direction in the 

Note that in these three acceleration examples the energy density distribution is flat around the origin 
r = 0 and therefore there is no singularity and, moreover, no cusp behavior (or weak singularity). 
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TABLE III: Three examples of the domain acceleration. 



t 

rD 

do 

rk 

nk 

hk 

n 

nt 

htb 

Qd 

1 

0.1 

1 

1 

0.6 

20 

10 

0.6 

20 

10 

- 0.01 

2 

0.1 

1.1 

10 ^ 

0.9 

40 

40 

0.9 

40 

10 

- 1.08 

3 

10 “® 

1 

lOio 

0.77 

100 

100 

0.92 

100 

50 

- 6.35 



Ld 

Ld 

Ld 

< 1 D 

1 

16.2 

1.62 

0.00174 

- 0.01 

2 

94.0 

7.63 

0.694 

- 1.08 

3 

8720 

117 

10.0 

- 6.35 


inhomogeneous region. How the evolution of the line elements affect the evolution of is 
not clear. Naively, it is a reasonable possibility that the domain acceleration stems from the 
acceleration of the radial line elements in the inhomogeneous region, that is, the existence 
of the accelerating line elements might be a necessary condition for the domain acceleration. 
Is it possible to have the domain acceleration without the acceleration of line elements? So 
far we do not have a no-go theorem prohibiting this possibility. However, we found no such 
example through our survey of the nine-dimensional parameter space. This might be an 
indication of the correlation between the domain acceleration and the acceleration of line 
elements. 

For a demonstration of the scales of the size and other quantities of this system, in Table 
IV we £x the only one unspecified unit by using the length unit: O.lMpc, IMpc, lOMpc 
and lOOMpc, and present the values of several dimensionful quantities (in addition to the 
dimensionless qo), respectively. For each example in Table HTTl when the size of the system 
increases by one order of magnitude, the time increases by one order and the energy density 
decreases by two orders of magnitude. From this table one can see that the time t and the 
energy density of the outer region cannot be simultaneously consistent with the situation of 
the present universe (i.e., fo ~ 10^°years and po ~ 10“^® g/cm“^), and therefore these three 
examples by themselves cannot describe the present universe. The search (based on the 
LTB solution) of the domain acceleration examples which are consistent with observational 
results is important and worthy of further investigations. 

To study the dependence of the deceleration parameter qjj on the nine parameters in 
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FIG. 5: The plots of the physical energy density p and the quantities, and which 

characterize the local acceleration/deceleration status in the radial and the angular direction, 
respectively, for the first example in Table Hill Note that when r = r/j = 1, dfy/g^ = —4.9 x 10“® 
and = —1.6 x 10“^ . 
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FIG. 6: The plots of the physical energy density p and the quantities, and d^^gee, which 

characterize the local acceleration/deceleration status in the radial and the angular direction, 
respectively, for the second example in Table Hill Note that when r = r£) = 1.1, df^/grr = —0.055 
and -0.43. 
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FIG. 7: The plots of the physical energy density p and the quantities, and df^gge, which 

characterize the local acceleration/deceleration status in the radial and the angular direction, 
respectively, for the third example in Table Hill Note that when r = vd = I, d^^/g^ = —0.23 and 
df^=-2M. 
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FIG. 8: Illustration of the dependence of the deceleration parameter qd on the parameters 
pQ,rk,nk, hk,rt,nt, htb), using the second example in Table Hill as a reference (denoted by 
the large dot). 
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TABLE IV: Corresponding to different length units, the values of several dimensionful quantities 
for the acceleration examples in Table IIIII are presented. 



<1D 

Length unit 

t (year) 

Lr (Mpc) 

p{r = 0) (g/cm^) 

p{r = ro) {g/cvD?) 

1 

-0.01 

0.1 Mpc 

3.26E4 

1.62 

3.15E-30 

4.81E-37 



1 Mpc 

3.26E5 

16.2 

3.15E-32 

4.81E-39 



10 Mpc 

3.26E6 

162 

3.15E-34 

4.81E-41 



100 Mpc 

3.26E7 

1620 

3.15E-36 

4.81E-43 

2 

-1.08 

0.1 Mpc 

3.26E4 

9.4 

2.11E-30 

1.79E-34 



1 Mpc 

3.26E5 

94 

2.11E-32 

1.79E-36 



10 Mpc 

3.26E6 

940 

2.11E-34 

1.79E-38 



100 Mpc 

3.26E7 

9400 

2.11E-36 

1.79E-40 

3 

-6.35 

0.1 Mpc 

3.26E-3 

872 

2.11E-16 

8.33E-36 



1 Mpc 

3.26E-2 

8720 

2.11E-18 

8.33E-38 



10 Mpc 

3.26E-1 

87200 

2.11E-20 

8.33E-40 



100 Mpc 

3.26 

872000 

2.11E-22 

8.33E-42 


Eq. fjdTl) , we use the second example in Table IIIII as a reference and tune one of these nine 
parameters at one time, while keeping the other eight unchanged (i.e., with the values in 
the second example in Table [In|) . The results are shown in Fig. |H1 These plots show that 
every parameter has signihcant influence on qd (at least within some range of the value of 
the parameter). We note that in the case of the inhomogeneous function tb{r) in Eq. fl38|) we 
have acceleration or deceleration in different situations with different values of parameters, 
and, moreover, the examples of the domain acceleration are not rare. Conversely, in the 
case of trivial tb{r) we found no domain acceleration without singularity. Accordingly, in 
the LTB solution the function tb{r) plays a key role in generating the domain acceleration. 

VI. SUMMARY, DISCUSSIONS AND OUTLOOK 

Against the common consensus [in Eq. ([T])] that normal matter always slows down the 
expansion of the universe, we have found and demonstrated examples of the line acceleration 
and the domain acceleration for a spherically symmetric dust fluid described by the LTB 
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solution without singularity. This discovery contradicts the common intuition about the 
interplay of gravity and the cosmic evolution. Furthermore, these examples have shown the 
strong correlation between acceleration and inhomogeneity. These results strongly support 
the suggestion that inhomogeneity can induce acceleration. 

In these acceleration examples the spherically symmetric dust fluid consists of three 
regions: two (roughly) homogeneous regions — the inner over-density region and the outer 
under-density region — and one transition/junction region, where the inhomogeneity locates, 
of two homogeneous regions. Note that in these examples the energy density distribution 
is flat around the origin r = 0 and therefore there is no singularity and, moreover, no cusp 
behavior (or weak singularity). 

For further understanding the acceleration, the quantity which characterizes 

the acceleration/deceleration status of an inhnitesimal radial line element, is one of the 
good quantities to study. Naively, the regions with positive/negative make posi¬ 

tive/negative contribution to acceleration. For the quantity d^^/g^ to be positive we find 
a necessary and sufficient condition: [a?/r)^r > 0. This condition tells us that a positive 
contribution to acceleration is made in the place where a(t, r) increases sufficiently fast with 
the radial coordinate r. In every acceleration example we found, there exists a region with 
positive d'fyfg^ that coincides with the inhomogeneous transition/junction region quite well. 
This result reveals the strong correlation between acceleration and inhomogeneity, thereby 
giving a strong support to the suggestion that inhomogeneity can induce accelerating ex¬ 
pansion. 

It is not clear how the line acceleration and domain acceleration are related to the ap¬ 
parent acceleration. A similar doubt (regarding the domain acceleration) was also raised by 
Kai et ah in Ref. 1^. The relation between the theoretical accelerations — line acceleration 
and domain acceleration — and the apparent acceleration could be model-dependent. For 
our universe with complicated energy distribution the relation may be far from simple. As a 
reminder, we note that the conclusion about the existence of the present cosmic acceleration 
indicated by observations is based on the homogeneous and isotropic FLRW cosmology. If 
one invokes inhomogeneous cosmology, the current observational results may not indicate 
the existence of the cosmic acceleration. As shown in Appendix [B] and Ref. 1^, in the cases 
where E{r) = 0 there is no acceleration (for both the line and the domain acceleration). 
In contrast, it has been shown that it is possible to £t supernova data in the LTB model 
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with trivial E{r) 2^. This gives an example of having apparent acceleration while having 
no line acceleration and no domain acceleration. One may treat a complicated universe as 
a large domain consisting of many different sub-domains, and from a statistical perspective 
it might be reasonable that the size of the large domain corresponds to the scale factor 
in the FRW metric. In many works the results from the analysis involving the quantity 
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60l |. Nevertheless, whether 


are compared with observations in this spirit 3y, [371, 1^, 
this is a good approximation is not clear yet. 

One may wonder whether the acceleration in our examples is real or fake, i.e., whether 
it truly corresponds to the space expansion. One example leading to this concern is to 
consider an accelerating system consisting of two decelerating regions. Since each region 
is decelerating one might expect that the domain acceleration in this case does not corre¬ 
spond to physically observable attributes (e.g. the luminosity distance-redshift relation) of 
an accelerating FLRW model 39j. Nevertheless, on the contrary, it has been shown that 
for a system consisting of decelerating regions it is possible to fit the supernova data or to 
generate apparent acceleration |2fil. l48[|. 

In addition, as already pointed out in Introduction, it needs much caution to connect 
separate regions or to put them together. Two particular important issues (that may be 
closely related to each other) are: 

(1) There may be singularity in the junction between separate regions. 

(2) The effect from the junction between separate regions may be signihcant and should be 
taken into account seriously. 


We have already had a lesson from the work by Nambu and Tanimoto in Ref. 1491, where 


they basically considered a system (of spherical symmetry) consisting of an inner and an 
outer decelerating FLRW region and, as a result, they found and presented examples of the 
domain acceleration. In this system there is singularity in the junction between these two 
regions, and in their treatment the contribution from the junction is ignored. As pointed 
out in the previous section (with details in Appendix O), it is inappropriate to ignore the 
junction. Actually the junction makes signihcant contribution to the volume of the system, 
although it looks like a “2D” surface in the coordinate space. After taking care of the 
singularity and taking into account the volume of the junction, we hnd that there exists no 


domain acceleration in the cases studied in Ref. 
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Another concern is about the dehnition of acceleration, as discussed in Sec. m The two 
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definitions employed in the present paper — line acceleration and domain acceleration, on 
which our acceleration examples are based — follow the scenario proposed in Sec. [TTl for a 
system of freely moving particles interacting with each other only through gravity. For the 
purpose of truly representing the space evolution status while avoiding the confusion from 
particle motion and fake frame acceleration, the length quantities invoked in this scenario 
to dehne acceleration are (1) the distance between two freely moving particles for the line 
acceleration and (2) the size of a spatial region with constant number of particles therein 
for the domain acceleration. Regarding the frame/gauge choice, we have emphasized the 
advantage of the synchronous gauge in which the above requirements for the length involved 
in the acceleration definition can be easily met. In addition, in this gauge, there exists a 
universal cosmic time that is the proper time of comoving observers, and for a region with 
its boundary hxed in the coordinate space the volume expansion rate coincides with the 
expansion rate of the physical proper volume. 

The counter-intuitive acceleration examples we found raise two issues worthy of further 
investigations: 

• How to understand these counter-intuitive examples? 

• Can inhomogeneities explain “cosmic acceleration”? 

Regarding the hrst issue, the common intuition about Eq. ([T]) may actually stem from 
Newtonian gravity that is a perturbed version of general relativity in the Newtonian limit 
with the Minkowski space-time as the background. Accordingly, this intuition may be 
valid only for considering the particle motion relative to a background space-time in a 
perturbative framework, but may be invalid for considering the evolution of space-time that 
is described by general relativity. For further understanding the cosmic evolution and other 
topics involving the general-relativity effects as the dominant effects, it will be helpful if one 
can build the intuition about general relativity, i.e., with which one can make a proper guess 
at the behavior of the space-time geometry for an arbitrarily given energy distribution. 

Regarding the second issue, the most important is whether the inhomogeneities of our 
universe can explain the supernova data. If yes, the next step is to see, according to obser¬ 
vational data, how the universe evolves, in particular, whether the cosmic acceleration exists 
or not, in this cosmological model taking the real inhomogeneities into account. So far the 
examples we found may be far away from the real situation of our universe. How to beneht 


from these mathematical examples in order to understand the present cosmic acceleration 
through the inhomogeneities of our universe is an important issue currently under our in¬ 
vestigations. No matter whether the comic acceleration can eventually be explained simply 
by inhomogeneity, these examples have shown that inhomogeneity may affect the cosmic 
evolution in a manner far beyond the usual naive intuition about the interplay of gravity 
and the space-time geometry. Thus, the effects of inhomogeneity on the cosmic evolution 
should be restudied carefully. These examples open a new perspective on the understanding 
of the evolution of our universe. 
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APPENDIX A: RESTRICTIONS ON THE LTB SOLUTION 

In this appendix, we will give several reasonable restrictions on the LTB solution and 
study the properties of the LTB solution under these restrictions. We consider three condi¬ 
tions: (1) no hole, (2) no singularity, and (3) non-negative and hnite energy density. For no 
hole at the center, the area of the spherical surface at r = tq should go to zero when tq goes 
to zero. Accordingly, Condition (1) requires 

R(t,r = 0) = 0. (Al) 

Condition (2) requires the smooth behavior of the functions involved in the LTB solution. 
In particular, we consider 

R,rit,r)^0. (A2) 

Condition (3) requires 

0 < p{t, r) < cx). (A3) 
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Note that the necessary and sufficient conditions of no shell crossing in a period of time are 


59l . The conditions discussed here are only necessary 


presented by Hellaby and Lake in Ref. 
conditions. 

Because is nonvanishing and continuous, R,r cannot change its sign, i.e., R,r is either 
always positive or always negative for all r. Because R(t, r = 0) = 0, the sign of R{t, r > 0) 


should be the same as that of R, 


i.e., 


RR,r >0 for r > 0 . 


(A4) 


This property will be used in Appendix O regarding the possibility of generating the domain 
acceleration with trivial tb{r). 

Considering Eq. ffT3|) . 

'r 

R 


2E(r) ^ 2M(r) 


R2 


R3 


for every term therein to be hnite, we obtain 


E{r = 0) = 0 , 
M(r = 0) = 0, 


(A5) 

(A6) 


because R{t,r = 0) = 0. From Eq. (ITTll . 


Eq. flASp requires 


p(t,r) = 


M'ir) 


2M'{r) 


> 0 . 


Because R,r cannot change its sign, M'{r) cannot change its sign either according to the 
above relation. Using the same reasoning for obtaining Eq. flA4p . from Eq. flA6l) we have 


M{r)M'{r) > 0. 


(A8) 


To sum up, so far we have obtained the following basic restrictions on the LTB solution: 
R{t,r = 0) = E{r = 0) = M(r = 0) = 0; R, R,r, M{r) and M'{r) have the same sign when 
nonvanishing. 

Examining Eqs. flT^ and fflTl) . we can see the feature that for a solution where r) < 0 
there is always another solution r) > 0, obeying 

Ri(f, r) = —R{t,r), Mi(r) = —M{r), Ei{r) = E{r). (A9) 
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Thus, without losing generality, we can set 


R{t, r) > 0, R^r > 0, M(r) > 0, M'(r) > 0. 


(AlO) 


With this choice, regarding the dependence on the radial coordinate r, i? is a monotonically 
increasing function and M is a non-decreasing function, while both vanish at the origin 
r = 0. 

In the following we will study -R(t, r) and show that R{t, r) < 0 and R{t, r = 0) = 0 
corresponding to the choice in Eq. flAlOD . The time derivative of Eq. flT^ multiplied by R^ 
gives a simple formula for R\ 

A/f ( 

(All) 


fl=-^hh<o. 




corresponding to the choice in Eq. (lAlOD . According to Eqs. flAlD and flA6p and applying 
L’Hospital’s rule, we have 


r = 0) = — 


M 


M' 


r=0 


2RR, 


r=0 


= --R.P 


(A12) 


r=0 


where Eq. ffTT|) has been applied in the last equality. As a result, for a hnite physical energy 
density p, 

ii(t,r = 0) = 0, (A13) 

because R(t,r = 0) = 0. 

APPENDIX B: NO ACCELERATION WITH E{r) = 0 = k{r) 


In this appendix we consider the trivial case, E{r) = 0 = k{r), and prove that there 
is no line acceleration and no domain acceleration in this case. (For the case of domain 


acceleration, see also Ref. 
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Line Acceleration 


In the case with trivial E{r) or k{r), Eqs. fl29|) and (l30|) give 

/•I’L rrL 

Lr{t)= / y/^rdr = / R,rdr = R{t,rL) , 
Jo Jo 

where Eq. (lAip . R{t,r = 0) = 0, has been used, and therefore 

RR 


Qr = 




R2 


(Bl) 


(B2) 


r=rL 
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From Eq. (fT^ . we have 


and 




2M(r) 

R 


R 


M{r) 


Substituting these two equations into Eq. flB2l) gives a simple result: 


Qr 


1 

2 ■ 


(B3) 

(B4) 


(B5) 


Thus, in the case with trivial E{r) or /c(r), no matter how we tune the other two functions 
Poir) and R{r), the deceleration parameter corresponding to the proper distance between 
the origin and any other point in space is always 1/2, a positive constant denoting a decel¬ 
eration as large as that in a homogeneous dust-dominated or pressureless-matter-dominated 
universe. 


Domain Acceleration 


When E{r) = 0 = /c(r), Eq. 

Vd = dvr 


gives 


'0 


R^R,r dr = td) , 


Lu = vI'^^(E^ ^ R(t,rD), 


and therefore 


Qd = 


LnL 


D-^D 


RR 


LI 




1 

2 ’ 


(B6) 

(B7) 

(B8) 


r=ro 


where Eqs. fIBSp and flB4p have been applied in the last equality. 

Thus, the deceleration parameters corresponding to the domain acceleration and the line 
acceleration are both 1/2, a positive constant (denoting deceleration) independent of the 
other two functions po{r) and This indicates that a non-trivial E{r) or k{r) function 

must play an essential role in the inhomogeneity-induced accelerating expansion based on 
the LTB solution. 
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APPENDIX C: DOMAIN ACCELERATION WITH TRIVIAL tb{r)l 


In Nambu and Tanimoto claimed that the example of the domain acceleration was 
found. However, we hnd that it is an improper example. In this appendix we will study the 
example in and discuss its problems. 


1. Nambu and Tanimoto’s Example 


In 


49l | the choice of the three arbitrary functions involved in the LTB solution are specihed 


as follows: 


Po(’") = Po = constant, (Cl) 

h{r) = 0, (C2) 

fc(r’) = ^ [26'(r - rfc) - 1], 0 < r < L, 0 < < L, (C3) 

where 9{r) is a step function. In the following study we will choose L = 1 without losing 
generality, that is. 


k{r) = 26{r — r^) — 1, 0 < r < 1, 0 < < 1. 


(C4) 


With the above choice the spherically symmetric dust fluid described by the LTB solution 
seems to consist of two regions described respectively by two different Robertson-Walker 
(RW) metrics: the inner region (0 < r < r^) described by a (spatially) open RW metric 
with the scale factor ai{t) and the outer region (r^ < r < 1) described by a closed RW metric 
with the scale factor 02 (t). The scale factors ai{t) and 02 (t) obey the following Friedmann 
equations, respectively: 


The initial condition 



ai{t = 0 ) = a2(t = 0 ) = 0 


(C5) 

(C6) 


(C7) 


are used to solve the above equations and to obtain the time evolution of the scale factors 
(oi and 02 ) and accordingly the time evolution of the domain {0 < r < L = 1}. 
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To obtain the deceleration parameter qr, of this domain, the following formulae for cal¬ 
culating the volume of the domain are used: 

(a -|- a,r r) dr 


Vnit) = dvr 


0 a/1 - k{r)r‘^ 

= dvr [cial{t) + C 2 al{t)] , 


(C8) 


where Ci and C 2 are constants dehned as follows: 

x'^dx x‘^dx 

Jo vT+~z? ’ ^ 

This volume calculation (i.e., the total volume being equal to the sum of the volumes of 
the inner region and the outer region) looks reasonable, but actually is incorrect, as to be 
discussed later. 

For demonstration, in the following we consider the special case with po = 3. The result 
about the time evolution of qc, obtained by using Eqs. flC5D - flC9D . is shown in Fig. [9l This 
hgure illustrates the existence of the domain acceleration at later times for several different 
values of r^. 




FIG. 9: Evolution of the deceleration parameter of a spherically symmetric domain, which 
consists of an inner open RW region and an outer closed RW region, for po = 3 and = 0.7,0.8,0.9. 

Note that with the initial condition in Eq. fICTIl we have ai{t) > a 2 {t) as t > 0. This con¬ 
tradicts the restriction in Eq. (lAdp around r = (i.e., around the junction of two regions), 

and accordingly indicates the existence of a singularity that stems from the increasing be¬ 
havior of k{r) around [but is irrelevant to the singular behavior of k{r) around r^]. This 
singularity can be avoided by exchanging the inner region and the outer region, making the 
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FIG. 10: Evolution of the deceleration parameter of a spherically symmetric domain, which 
consists of an inner closed RW region and an outer open RW region, for po = 3 and = 0.7,0.8,0.9. 


function k(r) decreasing instead of increasing, which is what we did about the choice of k(r) 
in our search for acceleration examples. 

Accordingly, with this change, k(r) becomes: 

k(r) =-2d(r - rk) + 1, 0 < r < 1, 0 < < 1. (CIO) 


In this case, the inner region {0 < r < is described by the closed RW metric with the 
scale factor ai(t) and the outer region {r^, < r < 1} by the open RW metric with 02 (t), 
satisfying a 2 (t) > ai(t) as t > 0. The result about qD(i) for the case with above k(r) and 
Po = 3 is demonstrated in Fig. [101 

As shown in Figs. [9]and[ini for both choices of k(r), there exists the domain acceleration 
(at later times). However, as to be shown in the following section, after smoothing the step 
function invoked in above k(r) in order to get rid of the other singularity, we hnd no domain 
acceleration, no matter how close to that in Eq. fIClOj) the function k(r) is. 


Problems and Revision 


In this section we use the following function k(r) to replace that involving a step function 
in the previous section. 


k(r) = 


2(r/rfc) 


nk 


+ 1 . 


(Cll) 


1 + (r/rk)^k 

When the power rik goes to inhnity, the above function k{r) approaches to that in Eq. (lCin|l 
involving a step function. If there is nothing wrong, we should obtain the same result about 
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goit) for the case invoking above k{r) with > cx) and the case invoking k{r) in Eq. fIClOD . 
For demonstration we choose t = 1, po = 3 and = 0.9, corresponding to the case where 
qd ~ —1 as shown in Fig. [10] for k{r) in Eq. fIClOD . In Fig. [TTjthe dependence of qd on 
for the case with k{r) in Eq. fICllD is illustrated. We can see that along with the increasing 
of the power the domain deceleration parameter always keeps far away from the value 
— 1 (that corresponds to signihcant acceleration), and smoothly approaches to a positive 
value around 0.16 (deceleration). 

0.3[^.: 

0.28 • : 

0.26- ■ : 

pO.24- . : 
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0.18- ’ ■ ■.. : 
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FIG. 11: The dependence of the deceleration parameter qo on the power Uk in the function k{r) 
in Eq. (jClip for t = I, po = 3 and = 0.9. 


Where does this discrepancy come from? The discrepancy stems from the doubtful cal¬ 
culation of the volume Vd in 49|] described in the previous section. In 49|] Nambu and 
Tanimoto ignored the volume of the junction (between two RW regions) at r = that is 
actually signihcantly nonzero (even though r = corresponds to a “2D surface” in the 
coordinate space) and should not be ignored. After taking back the volume at r = into 
account as a more appropriate treatment, we found no domain acceleration in the cases 


discussed in 


49|. 


In the following we study the volume of the junction region around r*, for the function 
k{r) in Eq. fICllD : 

rrk+e 


hj unction 47r 




— k{r)r‘^ 


where 


£ = 


2nk 


(C12) 


(C13) 


The dependence of Inunction, footai = Vd + Inunction and Inunction/Votai ou the power Uk is 
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illustrated in Fig. [121 This figure shows that the contribution to the total volume from the 
junction region remains significantly nonzero when Uk becomes large. 
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FIG. 12: The dependence of the volume of the junction and the total volume on the power in 
the function k{r) in Eq. (jCllI) . 


When Uk goes to infinity and accordingly e goes to zero, one can obtain a lower bound 
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of f^unction aS folloWS. 


rrk+e „ 2„2 


Uk 


OO ) unction 


Att 


> 


a^r^{a,r r)dr 

Jvk-e 
47ra?r? 

' a,r dr 


"v/l + 


J Vk € 


ATialrl{a2 - ai) 


V^TTr 


>0 if 02 > Cil , 


where £ = {2nk) ^ and 


(C14) 


ai(t) = lim a{t,rk — e), (CIS) 

rife—>oo 

a 2 (f) = lim a{t,rk + e). (C16) 

rife—>oo 

[Note that the hrst line of Eq. (lC14p is obtained by keeping only the singular part in the 
integrand.] Thns, even when the jnnction region goes to a “2D snrface” in the coordinate 
space, its volnme remains signihcantly nonzero if the scale factor of the onter open nnder- 
density region ( 02 ) is signihcantly larger than that of the inner closed over-density region 
(oi). As mentioned in the previons section, with the initial condition in Eq. flC7l) we do have 
02 > Oi for t > 0. 

Note that by changing the initial condition for oi(t) and 02 (f) [instead of nsing that in 
Eq. flC7[) ] Nambn and Tanimoto’s treatment conld become valid dnring some specihc time 
period when oi ~ 02 , which might be achieved throngh, for example, choosing non-trivial 
tb{r) (i.e., making the big bang or the expansion at different places begin at different times). 
This may be the reason why non-trivial tb{r) is a necessary ingredient in onr search for 
domain acceleration. 
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